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Abstract. We present a simple method for deriving the renormalization counterterms from 
the components of the energy-momentum tensor in curved space-time. This method allows 
full control over the finite parts of the counterterms and provides explicit expressions for each 
term separately. As an example, the method is used for the self-interacting scalar field in a 
Friedmann-Robertson- Walker metric in the adiabatic approximation, where we calculate the 
renormalized equation of motion for the field and the renormalized components of the energy- 
momentum tensor to fourth adiabatic order while including interactions to one-loop order. 
Within this formalism the trace anomaly, including contributions from interactions, is shown 
to have a simple derivation. We compare our results to those obtained by two standard 
methods, finding agreement with the Schwinger-DeWitt expansion but disagreement with 
adiabatic subtractions for interacting theories. 
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1 Introduction 

Quantum field theory in curved space-time has proven to be successful for studying the cos- 
mological implications of quantum fields, especially when interactions can be neglected [1, 2]. 
Other areas of physics have also benefited from the use of quantum fields in curved space [3] . 
Only recently have the effects of interactions made their way into such considerations. It is 
well-known that renormalization of interacting theories gives rise to important phenomena, 
such as the running of the couplings and in order to derive the quantum corrected equations 
motion, one must be able to perform consistent renormalization in curved space-time. The 
early work in this field concentrated on the consistent cancellation of divergences, without 
explicitly calculating the finite remainder of the counterterms [4-9]. However, if one is inter- 
ested in quantitative physical predictions, the finite parts of the counterterms, in particular 
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the renormalization scale, must be known in order to assign a proper physical interpretation 
for the constants of the theory. 

The various approaches to this problem can roughly be divided into two categories: 
Renormalizing at the level of the action and renormalizing at the level of the equations 
of motion. In the latter approach, the most popular method by far has been adiabatic 
subtraction [10-12] (see also [13]), having recently been applied in [14-16]. The former 
approach has also been found fruitful [17-20]. 

The ground work for renormalizing the quantum corrected action was laid in the seminal 
paper [21] and can be generalized to curved spaces [9]. Since the energy-momentum tensor 
is calculated as a variation with respect to g^ v , the effective action must be derived with a 
general metric. This is a highly involved task and there is currently quite a limited number 
of methods for deriving the effective action in curved space. Amongst the most popular is 
the gradient expansion, also known as the Schwinger-DeWitt expansion [22, 23], which was 
generalized in [24]. The gradient expansion works well for calculating the divergent contri- 
bution to the counterterms, but it is only reliable for the complete effective action when the 
fields are slowly varying, thus having limited applicability in many physical scenarios. An- 
other proposed method can be found in [25] (and references therein) but it is mathematically 
involved. 

The major benefit of adiabatic subtraction is that because working at the equation of 
motion level, one is free to constrain the metric, which in many scenarios offers a significant 
simplification over working with a general space-time. However, adiabatic subtraction does 
suffer from a disadvantage: If one does not introduce a procedure for also fixing the finite 
parts of the counterterms, the renormalization conditions for all the coupling constants is left 
implicit leaving the physical interpretation of these constants ambiguous. It has been shown 
that for the non-interacting scalar theory, adiabatic subtraction is equivalent to redefining 
the constants of the original action [26]. We are not aware of work examining the validity of 
adiabatic subtraction for interacting theories; whether the renormalization procedure can be 
reduced to the introduction of counterterms in the usual way. 

The aim of this work is to present a consistent and systematic renormalization procedure 
working at the level of the equations of motion, with explicit renormalization conditions for 
each constant and full control over the various finite parts of the counterterms. This method 
combines the benefits of adiabatic subtraction and the effective action approach. We will 
apply this method in the context of the self-interacting scalar field, i.e. a A</? 4 -theory, in a 
Friedmann-Robertson- Walker (FRW) space-time with zero spatial curvature, calculated to 
one-loop order. Explicit results will be derived for the energy-momentum tensor and the 
equation of motion for the (/j-field in the fourth order adiabatic vacuum. The regularization 
method will be dimensional regularization. The conformal anomaly, including contributions 
from interactions, is shown to have a simple derivation within this formalism. Finally, the 
results are compared to those obtained by the gradient expansion and adiabatic subtraction. 

The paper is organized as follows. In section 2 we discuss renormalization in curved 
space-time, with special attention being given to renormalization of the effective action and 
adiabatic subtraction. Section 3 contains the derivations for the one-loop equations of motion 
for the field and the metric in the adiabatic vacuum. In section 4 we present the explicit 
results for the equations of motion and the conformal anomaly. Comparison between our 
results and those derived via the gradient expansion and adiabatic subtraction is done in 
section 5 and we conclude in section 6. We will be using natural units where c = H = 1 and 
the (+,+,+) convention of [27] for the various tensors. Bare quantities are denoted with a 
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subscript " 0" , n — 1 dimensional spatial parts of vectors with a boldface letter and operators 
with a "hat", \ 



2 Obtaining renormalized equations of motion 

Quantum field theory in curved space-time is a theory of quantized fields on a curved clas- 
sical background [1, 2]. In the functional integral approach this means that the generating 
functional Z[J] has path integration over only the matter fields. For a theory with a single 
field ip and a classical action S[ip, g^ u ] this gives 



Z[J] 



Vip e 



(2.1) 



The quantum corrected or effective action denoted as Tlip^g^], can be derived by a 
Legendre transformation [28] 



TfagT] = J fx^g £ ef f&,gn = -HogZ[J] - j d 4 x^ J<p, 



(2.2) 



where it must be borne in mind that now (p represents the expectation value of the field, (<p) = 
(p. In general the effective action is divergent and renormalization must be implemented in 
order to have a meaningful result. For a renormalizable theory, the removal of the divergences 
is achieved by appropriately tuning the constants introduced by the original action Slip^g^]. 
After renormalization is done for T[ip, g^ u ], the finite equations of motion for the field and 
the metric can be derived via simple variation as 



5ip{x) 



0, 



5T[(p,g- 



5g» v (x) 



0. 



(2.3) 



One could equally well use the quantized action S[(p, g^ v \ to derive the operator equa- 
tions of motion 

/ ss[<p, 9 n \ = /ss[p,9n 







(2.4) 



\ 8<p(x) I ' \ 5g^{x) 

and renormalize the equations for the expectation values after the variation. Again for a 
renormalizable theory, we only need to adjust the constants of the action in order to render 
these equations finite. The first equation in (2.4) is the Schwinger-Dyson equation for the 
one-point function 1 [28]. 

In this paper our classical action will be 



S[<p,gr] = S m [<p,gr] + S g \gr\ 



W1 



2 ' ^ 



dntpd^ip + m 2 ip 2 + £Rtp 



d 4 x, 



-9 



A + aR + (3R Z + eiR a pR a/5 + e 2 R a ^sR 



2-^ 



(2.5) 
(2.6) 



We have included 0{R 2 ) tensors, which although not present in the classical Einstein-Hilbert 
action, are required in the quantized theory for consistently cancelling all the divergences in 
the energy- momentum tensor [2]. 



1 In this language, we will be truncating at one-loop, using bare vertices and propagators. 
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2.1 Renormalization via the effective action 



Given a regularized expression for the effective action r[<^,</" / ], renormalization can be ac- 
quired with relative ease: Each counterterm is determined by matching the effective action to 
the classical one, for each constant separately at some scale choices for the field, its derivative 
and the gravitational quantities. For an example, see [20]. In general this procedure must be 
used for all the constants of the classical theory and one must supply a scale for each degree 
of freedom. For the theory defined by (2.5 - 2.6) with a mean field that has a dependence 
only on time, the renormalization conditions are 



d 2 c, 



dip 2 
d 3 £ eff [if, g n 



1. 



d 2 £ eff [tp,g 



dip 2 dR 
d 2 C eff [if, g n 



dR 2 



(3, 



dip 2 



d 4 c, 



-m 



dR aP R a P 



A, 



dtp 4 
dC eff [ip,gn 



dC 



eff 



dR 

[<p,sT] 



dR a ^ 5 R a ^ 5 



-A, 



a. 



£2, 



(2.7) 

where tpi is a symbolic notation for all the matter and gravitational fields 2 and where all the 
scales Hi are arbitrary and can be chosen for each constant separately 3 . This procedure is 
simple to perform and it allows complete control over the finite parts of each counterterm 
separately. The effective action approach is however quite cumbersome for deriving the 
quantum corrected Einstein equation. If, for example, one is interested in the quantum 
corrections in a Friedmann-Robertson- Walker (FRW) space-time, the effective action must 
regardless be derived with a general metric and only after renormalization and variation with 
respect to g^ u as in (2.3) can the space-time be chosen to be homogeneous and isotropic. In 
general it is quite challenging to derive an expression for the effective action with an arbitrary 
metric. Usually one must settle for approximations such as the gradient expansion to be 
described in subsection 5.1. This has led to the construction of renormalization procedures 
operating at the level of the equations of motion. 

2.2 Adiabatic subtraction 

Adiabatic subtraction [10-12] is one of the most popular regularization (renormalization) 
methods used in curved space-time calculations, in particular for the energy-momentum 
tensor. The renormalized quantities are defined at the level of the equations of motion 
by subtracting an Ath order derivative approximation of the quantity of interest from the 
divergent expression. These quantities used in subtractions are said to be defined in Ath. order 
adiabatic vacua to be discussed in subsection 3.1. The order of the expansion A depends 
on the adiabatic order of divergences in the bare quantity. Thus in this procedure one only 
calculates one subtraction term that contains all the divergences needed for rendering a 
particular expression finite and counterterms such as Sm 2 or 5\ never explicitly enter the 
picture. For example, the renormalized variance of a field and the quantum part of the 
energy-momentum tensor defined via adiabatic subtraction read 

(P 2 ) = (v 2 )o + 5<f = (^ 2 ) - <0( A V|0^)| A=2 (2.8) 
(f&) = (T%)o + 6T„„ ee (fQ) - {0^\fQ\0^)\ A=d . (2.9) 



2 tpi = <p, tp2 = ip, ip3 = R, V>4 = -R a /3-R a ' 3 and ip s = R a/3 ^sR a 
3 Not all scale choices provide analytic solutions. See the discussion at the end of subsection (2.3). 
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The argument is that all the constants are taken to be renormalized since the divergences 
are taken care of by the adiabatic subtraction terms. 

This procedure has the advantage that it is generally easier to calculate the energy- 
momentum tensor starting directly from operator expressions than first deriving a renormal- 
ized expression for an effective action and then obtaining the energy-momentum tensor via 
variation. In particular, at the equation of motion level it is perfectly allowed to constrain 
the metric to be of special forms e.g. of the FRW-type. This method has its drawbacks, 
however. Because the renormalized expression is derived by a single subtraction, the renor- 
malization conditions for each constant including the renormalization scale, are not explicit. 
This leaves the finite part of each counterterm implicit and thus the physical interpretations 
for the renormalized constants become more difficult. The way to avoid this would be to add 
additional finite parts into the constants of the action and fixing them according to some 
chosen conditions, as was done in [29]. Since this fixing of finite parts ruins the calculational 
simplicity of adiabatic subtraction, it is rarely used. This is a rather unsatisfactory situation, 
since ultimately for a renormalizable theory one should be able to perform the entire renor- 
malization process by simply redefining the constants in the classical action. Furthermore, 
in [26] it was shown that for the case A = in (2.5), the adiabatic subtraction term in (2.9) 
can be reduced to a redefinition of the bare constants of the action, but it is a non-trivial 
issue whether this applies also to the interacting theory. This matter will be addressed in 
subsection 5.2. 

2.3 Consistent renormalization via the energy-momentum tensor 

In this section we set out to to find a renormalization procedure that allows full control over 
the finite parts of the counterterms while working at the level of the equations of motion, in 
particular the Einstein equation. This way one could combine the benefits of the effective 
action approach and adiabatic subtraction. 

Let us start by defining the energy-momentum tensor. For quantum matter and classical 
gravity denoted as S[tp, g^jo = S g \g^ v ]o + S m [<p, g^jo, 4 we get the Einstein equation via 
varying the bare action with respect to the metric 

2 6S[<p, 9 no =Q 2 6S g \gno = 2 £S m [y,<T]o =f (21Q) 

A crucial point is that although im>o has all the counterterms coming from the matter part 
via the bare coupling constants, we must also include counterterms from the gravitational 
side of (2.10) in order to remove all the divergences. Indeed, this is precisely the reason 
for introducing the higher order tensors in (2.6). We will denote the entire counterterm 
contributions as 

5T^^5T™-5T° U , (2.11) 

where 5TJ£ symbolises the counterterms from the matter action and 5T^ U the gravitational 
action. The renormalized energy-momentum tensor is then 



2_ 6S m [<p,g W 
~g Sg^ 



v = - * Tl + *v ( 2 - 12 ) 



4 In our notation the subscript "0" denotes a quantity where all the coupling constants are bare constants. 
For renormalized constants, we simply drop the subscript. A bare constant Co can be split into a (finite) 
renormalized part and a counterterm as Co = c + 5c 



- 5 - 



For the gravitational part in (2.6) with the help of appendix C we have 



2 5S g [g 



-9 Sg^ 



-g^A + 2a GV + 2/3 + 2e lj0 {2) H^ U + 2e 2fi H l 



flUl 



(2.13) 



so in addition to the standard counterterms from the matter part of the action, such as 5m 2 
and 5X, we have the gravity counterterms 



STjtv = -g^SA + 25aG^ + 25(3 Wfl^ + 28e x ^H^ + 25e 2 H l 



(2.14) 



We would now like to find a generalization of the equations (2.7) for the energy- 
momentum tensor, so that we can determine each counterterm in (2.12) separately. For 
this we can use equations very similar to those in (2.7). The only difference being that we 
use matching between the complete expectation value of the energy-momentum tensor to the 
classical one for determining the counterterms. We no longer need to keep a general metric, 
since the variation with respect g^ is already done. Therefore we can choose constraints 
for the metric, for example that we are in a homogeneous and isotropic space, and use this 
accordingly in the renormalization conditions. 

Now we will derive the renormalization equations for the coupling constants in (2.5) 
and (2.6) in a space-time with a FRW-type metric g flL ,dx^dx u = —dt 2 + a(t)dx 2 with a non- 
zero mass and a possible time dependence in the field <p. Since we will employ dimensional 
regularization in section 4, these formulae are calculated in n dimensions. Using appendix C 
we straightforwardly derive the result for the classical energy-momentum tensor 



rpC 



1 



in 



1 



1 



which gives the renormalization conditions 



00/ 



dip 2 
d 3 (f 00 ) 



d 2 (T, 



oo 



dtpdtpd(d/a) 
5 3 (f 00 ) 



d(h/a) 2 d(d/a) 



2£(n - 1) 



dip 2 
(Too) 



d 2 (T, 



oo 



d(d/a) 



m 
0, 

o. 



if 2 + 2(n — 1)— <ptp 
a 



(2.15) 



<9 4 (T C 



00/ 



A. 



d 2 (T, 



oo 



d(a/a) 2 
d 4 (f 



'-oo 



d{a/a) 



0. (2.16) 



Here again ip symbolizes the expectation value of the field just like in (2.7) and the scale 
choices for each constant can be done independently. The first four conditions in the above 
fix the counterterms coming from the matter part of the action and the rest do the same 
for the counterterms from the gravity part. It is obvious that there are many other possible 
choices for renormalization equations. For example 5£ can also be determined from 



00/ 



dp 2 d(d/a) 



4£(n 



1 



(2.17) 



Also, we could have chosen any other component of T^ v for determining the renormaliza- 
tion constants. Had we chosen to use different conditions to those in (2.16), the resulting 
counterterms potentially would have had different finite parts. 
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A few comments are now in order. First, since renormalization was achieved by simply 
including counterterms in the constants of the original action, covariant conservation is au- 
tomatically satisfied 5 . If the divergences are removed by introducing a subtraction term by 
hand, covariant conservation must be separately checked. Second, a crucial requirement for 
the equations (2.16) to work is that all the equations are analytic at the chosen scales /ij. As 
was discussed in [21] in the </? 4 theory there is an infrared singularity in the massless limit. 
To bypass this issue at the renormalization stage, one may choose non-zero renormalization 
scales fa. However, changing renormalization scales changes the definitions of the constants 
and ultimately changes the parameter space where the perturbative expansion is valid. An 
example of using a non-zero renormalization scale can be found in subsection 4.2. To cure 
the IR problem altogether requires further resummations [30]. 

3 Field equations and the energy-momentum tensor 
3.1 Choosing the vacuum 

For our calculations we will use the adiabatic vacuum [1] and we must perform our calculation 
in n dimensions since dimensional regularization is to be used in section 4. With a field 
equation of motion of the form 



U + M 2 {t) 



= 0, (3.1) 



where M(t) is some arbitrary and possibly time-dependent mass parameter, we can solve it 
in terms of mode functions by using an ansatz 

= [ d n ~ l k [a k u k + 44] , u k = — 1 =/t k (t)e^ x , 
J 1 J v / 2(2vr) ra - 1 a n - 1 

k k [t) = -$=e- i ?' Wd * (3.2) 

Vw 

and assuming g^dx^dx 11 = —dt 2 + a(t)dx 2 , we can write W as an adiabatic expansion 

a M a 2 M 2 a M dM . , 

W = c + ci- + c 2 — + c 3 — + c 4 — 2 + c 5 ~ + c 6 — + c 7 — + • • • , (3.3 
a M a A M z a M am 

with Ci being functions of M and a. This approximation can be trusted when ip and a are 

slowly varying. The solution for will naturally also be an expansion in the number of time 

derivatives and the Ath order approximate solution will include all terms with an A number 

of derivatives and it will be denoted as . Similarly, the vacuum it defines is written as 

|0^). These approximate modes can be used to define an exact solution to the equation of 

motion (3.1) through the relation 

n k = a k (t)^ A) +/3 k (t)4 A) *, (3.4) 

where a k and /3 k must be constant in t to order A, but they may have dependence in k. We 
can then choose to fix the exact mode ii k at some point t = to to be the ^4th order positive 
solution: 

a k (t ) = 1 + C(A+1) /3 k (£ ) = + 0(^ + 1) (3.5) 



5 This can be shown by operating with V M on the gravitational side of (2.10) and using the Bianchi identities 
and commutator formulae. 

6 This means that a 2 and a are of the same adiabatic order, for example. 
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This defines the Ath. order adiabatic vacuum. Clearly this procedure has some ambiguity, 
since we can choose an arbitrary point to as our condition to fix the exact mode resulting in 
different ^4th order vacua u k . However, the exact modes differ only in the ^4 + 1 adiabatic 
order, so if we choose to neglect 0(A + 1) we can write the quantized field operator to O(A) 
as 

.<^At..CA).l (3 .6) 



<p(x) = J dk n - x [a k u£ } + , 
with the standard commutation relations 

[flk, Ok'] = [a k , a k /] = 0, [o k , a k ,] = <5(k - k'), 



(3.7) 



and remember that in principle it is the exact mode that becomes quantized and defines 
the vacuum. Troughout this paper we will work in the adiabatic approximation, where we 
include only two or four time derivatives (A = 2 or 4). 

3.2 Field equation to one-loop order 

From the action in (2.5) the equation of motion for the operator ip reads 



+ ^ = 0. 



(3.8) 



Since we are interested in the one-loop corrections, we can shift the field operator as <p — > 
ip + 4>, where (<p) = (p as in section 2.1, and expand (2.5) around = giving to quadratic 
order 7 



d n xyf^g 4> 



4 



□ + ml + CqR + y¥> 2 



4! 

+ 



(3.9) 



In a similar fashion we can also split (3.8) into two coupled equations, one for the expectation 
value of the original mean field and one for the shifted operator 



□ + ml + Co-R 



3! 



a + ml + toR + Y 1 ? 2 



. A 3 , A ,- 2 



2 
0. 







(3.10) 
(3.11) 



Note that we also used the fact that since (j) is Gaussian at one loop order, all correlators 
with an odd number of fields vanish. Again, as indicated by the subscripts "0", we need 
to split our bare coupling into a renormalized part and a counterterm. Since we include 
only the first non-trivial correction and since the counterterms are the result of the one-loop 
quantum correction, all terms such as 5\((j) 2 } are effectively two-loop terms and beyond our 



7 In the one-loop approximation, the terms linear in cf> can be discarded. This is because we can always 
choose a renormalization condition for a linear term in the classical action that cancels this contribution [28] . 
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approximation 8 . Inserting the counterterms into equations (3.10) and (3.11) we get 



□ + m 2 + 5m 2 + (f + 5£)R 



A + SX o A . 79. 



□ + m 2 ^ 2 



0. 



(3.12) 
(3.13) 



In the above we used the fact that in a one-loop calculation there is no need for a wave 
function counterterm multiplying the d'Alembertian operator 9 . 
The equation for the operator (3.11) <fi 



d t d t + (n - l)-dt - aT 2 did 1 + m 2 + ^ip 2 + £R 
a 



0, 



(3.14) 



is precisely of the form (3.1) and hence we can use the ansatz from (3.2), which after some 
algebra gives an equation for W 



W 2 = a(t) + 



3W 2 W 



AW 2 W 



where 

a(t) = k 2 /a 2 + m 2 + ^ 2 +- 
2 a 



_(n-l)(4£-l) 



a 

+ 1 - 

a 



-(n-l)(3-re + 4(re-2)£) 



(3.15) 



(3.16) 



and the re dimensional scalar curvature in FRW was found from (C.10). W 2 can now be 
solved from (3.15) iteratively as 

W 2 = (W 2 )^ + (W 2 p + (W 2 )^ + (W 2 )^ + . . . , (3.17) 

where, as explained in section 3.1, the expansion parameter is the number of derivatives. We 
will need W only up to four derivatives and the complete results are relegated to appendix 
A. 

3.3 Energy-momentum tensor to one-loop order 

From (3.9) we get the re-dimensional unrenormalized energy-momentum tensor to one-loop 
order by straightforward variation: 

-2 SS m [ipJ,g^] 



'-9 



Sg 



9tLV 

2 

9jiv 
2 



4 



4! 



+ d^ipdvtp + £ [Gpv - V^V„ + g^D] ip 2 
+ d^dj + Co [Gjju, - V M V^ + g^D] 4> 2 



rpC 1 n~>Q 



(3.18) 



s This amounts to counting loops at the level of the effective action, rather than at the level of the equations 
of motion. The discarded term would be included if the the two- loop "figure-8" diagram was included at the 
level of the action. 

9 Even if there is no divergence proportional to the kinetic term, one might still need a counterterm in order 
to give a precise physical interpretation for this quantity. For our renormalization scale choices this will be of 
no relevance. 
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That the energy-momentum tensor nicely splits into a classical and a quantum part is slightly 
misleading since only the complete has the right conservation properties. Including the 
counterterms from the gravity side and separating T^ v into renormalized and counterterm 
parts, we have for the renormalized E-M tensor 

T^^TC+fQ+ST^, (3.19) 

where the counterterm contributions for our one-loop approximation are 

+ 8£ [Gfiu - V^V U + g^U] p 2 - 8T*,, (3.20) 

where 5T^ U was defined in (2.14). 

In the four dimensional conformal limit we have m = and £ = 1/6. For our n 
dimensional expression we then have 

m 2 = and £o = ""^ + (3.21) 
4(n — 1) 

Inserting these into T^ v and (T^), we get the relations 

^% C ,,o = ^^^0 2 ) + (n-4)^^ +( 5e(n-ip^, sT {!$,) = -<f±(p), (3.22) 
giving for the full trace 

+ 5£(n-l)Op 2 -(5TV, (3-23) 

which course will result in the famous conformal anomaly first seen in [31]. The standard 
purely gravitational contribution is from the last term of (3.23) and the first three give the 
one-loop A-dependent terms, which are studied in more detail in [32, 33]. The equations in 
this section also provide a very non-trivial check for the solutions to be obtained in section 
4. 



ST, 



2 



X 2 2 , n^X 4 

dm p + 2 



4! 



(T/) = (n - 4) 



Xip 2 



4 Explicit results 

Next we will proceed to use the methods described in previous sections for the theory de- 
fined by (2.5) and (2.6) in the four dimensional case and present the result for the equation of 
motion for the field tp and the expression for the energy-momentum tensor calculated in the 
adiabatic vacuum as defined in subsection 3.1. A point worth emphasizing is that because all 
the calculations are done in n dimensions, everything is dimensionally regularized and hence, 
in practice, finite. For consistency it is important that one starts from an n-dimensional 
Lagrangian and calculates all the required quantities in arbitrary dimensions (as we have 
done), as opposed to working down from a 4-dimensional theory and only regularizing ex- 
plicitly divergent equations. We end the section by calculating the result for the anomalous 
four dimensional trace. We will notate the classical parts and the various adiabatic orders of 
the quantum contributions as superscripts 

(f^) = T° v + (f Q )(°) + (f^ + (f^ 
/ 6S f:^° \ =e° + £<°> + £ W + £ W = 0. (4.1) 

\ <w) / 
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4.1 Results to fourth adiabatic order 

Using the expressions from appendix C and subsection 3.1, we can write the non-trivial 
components of the quantum energy-momentum tensor in the adiabatic vacuum as 
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where M 2 = m 2 + \p 2 ■ By inserting (A.l - A. 3) into (4.2) gives the explicit result for "00" 
components of the quantum part of the energy-momentum tensor, which can be used to 
perform renormalization as explained in subsection 2.3 with the scale choices 
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where the condition a = 1 states that all derivatives of a are zero. By using (B.l - B.7) in the 
expression for the full T^ u in (3.19) calculated with (A.l - A. 3) and using standard formulae 
for n-dimensional integrals [28], for n = 4 we get 
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The same counterterms renormalize the "ii" component of the energy- momentum tensor, 
with T~ given by (4.3), and the equation of motion for the field tp (3.12). Similar steps that 
led to (4.5) give 
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and for the field equation of motion (3.10) 
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Covariant conservation was explicitly checked. We also see that at order and 2, the deriva- 
tives of a nicely organize themselves as powers of H = d/a and scalar curvature R. In general, 
it must be possible to write the result in terms of covariant tensors to the extent that the 
gradient truncation includes all derivatives of a given order. 

4.2 The conformal anomaly 

As was pointed out at the end of subsection 2.3, in the massless limit in a one-loop calculation 
we can formally work our way around the infrared singularity at m = by renormalizing A 



at a non-zero scale for the field <p. So in the conformal limit £ 
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and m = 0, we can 



use the counterterms defined in (B.l - B.7) with the exception of 5X, which we define as 
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which coincides with [34]. We point out that the conformal anomaly can also be derived 
from the renormalized four-dimensional components of T^ u . This calculation can be found 
in appendix D. 



5 Comparison with other methods 

5.1 Comparison with the derivative expansion 

Since we have used an adiabatic expansion for the modes in obtaining the results (4.6 - 4.10) 
it is illuminating to compare the results of the previous section to the equations of motion 
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derived via the gradient or Schwinger-DeWitt expansion [22, 23], which is also in essence 
adiabatic. It is believed that results derived via the Schwinger-DeWitt expansion for the 
effective action and the adiabatic expansion for the modes in fact give identical results. For 
a more detailed exposition of this method we refer the reader to [20]. 

We start by writing the standard loop expansion for the effective action 



r>, 5 H = J d n x^- g c eff = r«V<T] + r (1 V, 5 H + 
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where for (2.5) and (2.6) we have 

r(°V, 9 H = s m [<p )S r] Q + s B n, rW^H = ~Tr log g-\x, x '), (5.2) 

and G(x, x') is the Feynman propagator satisfying 
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The dimensionally regularized derivative expansion for the one-loop contribution reads 
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and the scale \x was introduced to maintain proper dimensions for the n-dimensional action. 
We will calculate the result only up to two derivatives (i.e. to second adiabatic order) and for 
this we need the Schwinger-DeWitt coefficients from ao to that can be found in appendix 
E. The result is 
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and discarding terms of 0(R 2 ) we can derive the counterterms, which can be found in ap- 
pendix E. Inserting these into (5.6) gives the effective Lagrangian 



C eff [<p,!r]=A + aR-± 



[9AV + ^ {2m 2 + m ~ 1)R) + 4^(6? - 1)R] + 
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where the last term symbolises the neglected fourth adiabatic order. From this expression 
the equations of motion can be derived by variation 
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If we choose = diag(—l, a 2 , a 2 , a 2 ) the results are precisely those in (4.5), (4.7) and (4.9) 
so to this order the Schwinger-DeWitt expansion for the action gives identical results to the 
adiabatic expansion of modes. It is also plausible that this equality holds as well for higher 
order results, but of course the proof of this assumption requires higher order calculations. 

5.2 Comparison with adiabatic subtraction 

In this subsection we will focus on the time components of the energy-momentum tensor. 
Let us first write the finite parts of the subtraction term 5T^ U from (3.20) by using the 
counterterms (B.l - B.7) 
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where the dots indicate contributions that vanish in the four dimensional limit. We have 
introduced an arbitrary scale \i for dimensional purposes and used it to absorb a factor of 4ir 
and the Euler gamma constant from the n-dimensional integrals. We explicitly see that in 
the one-loop approximation the finite part of the subtraction term STqq is formed solely from 
the terms in the classical action (2.5 - 2.6) i.e. it is a polynomial in the degrees of freedom, 
(p, a/a etc., as it is apparent from the derivation in subsection (3.3). 

In a similar fashion we can write explicitly the finite logarithmic parts for the energy- 
momentum tensor calculated in the fourth order adiabatic vacuum 
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All the finite terms not included in (5.12) are proportional to a finite inverse power of M 2 
and vanish at A = 0. In adiabatic subtraction —(Tq ) is precisely the term one uses to render 
the bare energy-momentum tensor finite as in (2.9). Comparing (5.11) with (5.12) we see 
that there is a discrepancy between the two expressions. In particular in (5.12) there are 
terms that are not polynomials in ip and hence cannot be obtained from the classical action 
by redefining its constants, which can be seen from the M 2 = m 2 + Ay? 2 /2 dependence of the 
logarithm. When one sets A = we have M — > m and the two expressions coincide verifying 
[26]. Hence we conclude that adiabatic subtraction coincides with the method described in 
subsection (2.3) only at the non-interacting limit, taking all renormalization scales to zero 
for the gravitational counterterms. In fact this result could have been easily deduced from 
the expressions in (4.6 - 4.10): We chose the fourth adiabatic order vacuum for the bare 
energy-momentum tensor, which means that in the context of adiabatic subtraction, the 
bare expressions are equal to the subtraction terms in (2.9) and hence adiabatic subtraction 
gives identically zero results for all the quantum contributions. This is the case for (4.6 - 
4.10) only when A = 0. 

6 Summary and Conclusion 

We have shown how to derive the renormalization counterterms consistently from the com- 
ponents of the energy-momentum tensor in curved space-time. This procedure also allows 
complete control over the finite parts of the counterterms. When working at the equation of 
motion level one is free to constrain the metric to be of the desired form. We then used our 
method for calculating, in the one-loop approximation, the renormalized equation of motion 
for the field and the components of the energy momentum-tensor to fourth adiabatic order in 
the adiabatic vacuum. We also calculated the result for the anomalous trace, in the one-loop 
approximation and finally compared our results to those obtained via two standard methods: 
The gradient expansion and adiabatic subtraction. 

We find that the method proposed here has advantages over both the gradient expan- 
sion and adiabatic subtraction for the reasons just described: Being able to constrain the 
metric allows us to specialize and solve for a wider range of models and explicit knowledge 
of the finite parts of the counterterms is important for the physical interpretation of the 
coupling constants. The steps described in subsection 2.3 are completely general and should 
be applicable to various different vacuum states and different theories. We chose to use our 
method for the adiabatic vacuum, since this provided fruitful comparisons to the gradient 
expansion and adiabatic subtraction. We were able to verify that our method gives identical 
results to the Schwinger-DeWitt expansion to second adiabatic order, which likely indicates 
an equivalence in higher orders as well. 

The comparison with adiabatic subtraction revealed that when interactions are included, 
adiabatic subtraction introduces counterterms that cannot be obtained from redefining the 
coupling constants of the theory alone. Our procedure has no such problems. For the non- 
interacting case this discrepancy is removed. 

When using this formalism in situations where the energy-momentum tensor involves 
convolutions over the various degrees of freedom, such as higher order loop expansions, ad- 
ditional complications may arise since the energy-momentum tensor is no-longer a simple 
function of the degrees of freedom, and a mode function anzats for the field may be inap- 
plicable. In addition, non-analyticities in the renormalization scale choices may pose further 
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obstacles. Nevertheless we believe that this method will be usable in various applications of 
curved space field theory due to its generality and mathematical straightforwardness. 
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A Solutions for the modes up to four derivatives 

From (3.15) we find the zeroth order solution of W 2 in (3.17) to be 
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where M 2 = m 2 + A/2^ 2 , w 2 = k 2 /a 2 + M 2 and 7 = M 2 /u 2 . 
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B Counterterms for the adiabatic energy-momentum tensor 

Using the renormalization conditions from (4.4) for the quantum energy-momentum tensor 
components (4.2) and (4.3) calculated with the fourth order adiabatic modes, we get the 
counterterms 
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where ui 2 = k 2 + m 2 . The surprising thing is that 5e2 is not needed: when 5j3 and <5e2 are 
inserted into 5Tf LV the coefficient of 5e\ vanishes identically. 

C Geometric tensors in n dimensional FRW spaces 



Standard variational calculus gives the following geometric tensors 
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and 
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The n dimensional formulae can be derived by noting that irrespective of the dimensionality, 
the non-zero components of the Riemann tensor in conformal time, 
g^udx^dx" = a(r])(—dr) 2 + dx 2 ), are 
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Using the above expressions we can now write our tensors in Minkowski time. The second 
order tensors 
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and the fourth order tensors 
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D Conformal anomaly from the renormalized T, 
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If we repeat the calculation of section 4, but with renormalizing the A at a non-zero scale as 
in (4.11) and set £ = 1/6 and m = we get 
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and 
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Denoting the left hand side of (D.3) as £, it can now be used to rewrite the trace of T„v in 
various equivalent forms. If we choose the following expression for the trace 
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we get exactly the from in (4.12). 



E Counterterms and Schwinger-DeWitt coefficients for the effective action 

The first few ao terms in the gradient expansion (5.4) can be found in e.g. [35] 
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where the dots indicate terms with an adiabatic order higher than 2. 

From the renormalization conditions in (5.8) we get the following results: 
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